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ABSTRACT 


A  technique  called  Nearfleld  Acoustical  Holography  (NAH) 
has  been  developed  at  the  Pennsylvania  State  University. 
The  technique,  a  significant  improvement  over  conventional 
holography,  is  not  subject  to  the  wavelength  resolution 
limitation  associated  with  conventional  techniques.  An 
exact  formulation  of  the  Greene’s  function  propagator  allows 
reconstructions  of  the  entire  pressure,  particle  velocity, 
and  vector  intensity  fields.  The  vector  intensity  may  be 
plotted  to  show  the  flow  of  acoustic  energy  around  and  away 
from  a  source.  On  some  sources,  these  mappings  reveal  areas 
of  non-radiating,  circulating  energy  flow.  These  areas  may 
have  large  pressure  or  velocity  amplitudes  and  negative 
intensity  regions  that  previously  led  to  misconceptions 
about  the  production  of  acoustic  energy. 

This  work  presents  details  of  the  development  of  the  MAH 
system.  Included  are  a  review  of  the  theory,  details  of  the 
experimental  apparatus  design,  accuracy  verification,  and 
the  use  of  the  system  to  study  a  relatively  simple  source  (a 
homogeneous  rectangular  plate)  and  a  complex  source  (an  in¬ 
tact,  string  excited  guitar).  The  guitar  research  revealed 
some  interesting  results  such  as:  the  rose  is  an  important 
source  of  low  frequency  energy,  the  motion  of  the  back  plate 
is  very  important,  and  similar  plate  modes  occur  at 
different  frequencies  due  to  a  complex  relationship  among 
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CHAPTER  ONE 

INTRODUCTION 

I.l  Statement  of  the  Problem 

Mechanical  vibrators  and  their  radiated  sound  fields  have 
been  of  considerable  Interest  In  the  past.  Yet  how  the 
motion  of  the  vibrator  couples  to  the  medium  to  produce 
sound  and  how  this  sound  energy  flows  away  from  the  vibrator 
to  the  farfleld  remains  Inadequately  understood. 
Mathematical  predictions  of  the  sound  field  produced  by  a 
vibrator  are  quite  Involved  and  only  a  few  Idealized  sources 
may  be  studied  In  this  manner.  The  mathematics  necessary 
for  studying  a  complex  vibrator  exceed  the  realm  of 
practical  solution.  Consequently,  experimental  methods  have 
been  pursued  to  study  complex  vibrators.  Many  existing 
measurement  techniques  accurately  record  the  properties  that 
they  were  designed  to  measure,  but  the  Information  provided 
to  the  researcher  does  not  reveal  the  total  radiation 
picture.  For  obvious  academic  and  practical  reasons.  It 
would  be  advantageous  to  develop  better  experimental 
techniques  and  expand  our  understanding  of  the  relationship 
between  a  complex  vibrator  and  Its  radiated  sound  field. 

The  purpose  of  this  work  Is  to  present  a  new  technique 
called  Nearfleld  Acoustical  Holography  that  has  many 
distinct  advantages  over  other  techniques  for  sound 
radiation  measurement.  This  >  tesls  '  .1  cover  the  design 

and  construction  of  the  experimental  apparatus,  accuracy 
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checks  on  the  system.  and  the  use  of  the  system  with  a 
simple  and  a  complex  vibrator. 

Nearfleld  acoustical  holography  Is  similar  to 
conventional  acoustical  holography  but  an  Important 
distinction  allows  much  more  Information  to  be  obtained  from 
the  measurement.  This  will  be  discussed  In  detail  later. 
One  Important  result  from  this  technique  Is  an  energy  flow 
map  that  shows  the  actual  magnitude  and  direction  of  the 
acoustic  energy  flow  field.  These  flow  maps  Indicate  which 
areas  on  a  vibrator  are  producing  and  absorbing  energy  and 
also  reveal  regions  of  circulating  energy  flow. 

Before  describing  the  differences  between  nearfleld  and 
conventional  acoustical  holography,  a  discussion  of  the 
basic  problem  Is  In  order.  One  would  like  to  be  able  to 
correlate  the  mechanical  properties  of  a  complex  vibrator 
with  the  radiated  sound  field  It  produces.  This  Is  not  an 
easy  connection  to  make.  Exact  mathematical  solutions  that 
predict  the  radiated  sound  field  are  available  for  only  a 
few  simple  sources.  If  the  source  Is  not  homogeneous,  as  In 
the  case  of  a  plate  with  stiffening  rib,  then  the  precise 
calculation  of  the  sound  radiation  becomes  extremely 
difficult.  If  not  Impossible.  The  next  best  solution  may  be 
to  experimentally  study  the  source  using  conventional 
acoustic  measuring  techniques.  But  these  measurements  may 
be  misleading.  Take  as  an  example  two  plates  vibrating  In 


some  natural  mode  with  one  of  the  plates  above  coincidence 
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and  the  other  below  coincidence.  Measuring  both  plates  with 
an  accelerometer  may  reveal  similar  surface  velocities  In 
both  plates,  but  the  plate  operating  below  coincidence  will 
not  radiate  much  energy  Into  the  farfleld.  The  plate 
operating  above  coincidence,  though  It  has  roughly  the  same 
velocity,  radiates  a  larger  amount  of  energy  Into  the 
farfleld.  Thus  a  knowledge  of  the  surface  displacement  or 
particle  velocity  around  the  vibrator  does  not  necessarily 
determine  how  sound  energy  Is  transmitted  to  the  sound 
field.  In  general,  the  problem  is  not  so  severe  at  higher 
frequencies,  above  coincidence,  where  large  surface 
displacements  most  likely  produce  the  radiated  energy.  The 
difficulty  occurs  at  lower  frequencies  where  the  radiated 
wavelength  Is  larger  than  the  source  dimensions.  At  lower 
frequencies  large  displacements  or  pressure  amplitudes  may 
actually  correspond  to  sinks  of  acoustic  energy. 

To  accurately  determine  how  a  vibrator  radiates  sound  and 
how  this  energy  flows  Into  the  farfleld,  a  knowledge  of  the 
acoustic  Intensity  vector  field  Is  required.  This  Is 
defined  as  a  product  of  the  In-phase  components  of  the 
pressure  and  the  particle  velocity,  l.e. 

Uh  -  1/2  p(r)  vCO  cos  0(0  (1) 

where  9  Is  the  phase  difference  between  the  pressure  and  the 
velocity.  Since  this  Is  a  vector  It  gives  the  rate  and 


direction  of  the  acoustic  energy  flow  at  any  point  in  space. 
A  large  component  of  this  quantity  normal  to  the  vibrator 
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corresponds  to  an  energy-producing  area  on  the  vibrator. 
Because  of  this.  Interest  In  measuring  t(r)  has  recently 
grown.  A  two-mlcrophone  technique  for  measuring  Intensity 
Is  available,  but  It  Is  limited  In  that  It  records  only  one 
component  of  the  Intensity  at  one  point.  This  limited  data 
can  lead  to  confusion  In  defining  the  energy-producing 
regions  that  radiate  to  the  farfleld  because  of  the 
existence  of  non-radlat Ing ,  circulating  Intensity  fields. 
Using  the  two-mlcrophone  technique  to  map  a  sufficiently 
large  energy  flow  field  would  be  quite  time  consuming. 

The  following  Is  a  list  of  some  of  the  Important  features 
of  nearfleld  acoustical  holography^: 

1.  One  non-contact  measurement  allows  analysis  of  the 
entire  sound  field  above  the  vibrator. 

2.  Data  taking  Is  accomplished  quickly,  requiring  roughly 
15  minutes  from  beginning  of  the  data  aqulsltlon  to 
analyzing  the  results.  This  allows  the  researcher  more 
time  to  study  the  results  rather  than  record  data. 

3.  The  measurement  area  Is  large  so  that  sources  can  be 

•) 

pinpointed  within  roughly  a  10  m~  area. 

4.  The  measurement  area  subtends  a  large  solid  angle  from 
the  source  and  therefore  can  measure  mul t 1-d 1 rec t lonal 
sources  without  missing  Information. 

5.  High  spatial  resolution  permits  localization  of 
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sources  to  within  about  10  ca. 

6.  Output  can  be  computer  graphic  displays  of  the: 

a.  sound  pressure  field  from  the  source  farfleld, 

b.  particle  velocity  field  from  source  to  farfleld, 

c.  modal  structure  of  the  vibrating  surface  (obtained 
from  the  normal  velocity  at  the  source  surface) , 

d.  vector  Intensity  field  -  this  locates  sources  of 
sound  energy  and  maps  the  energy  flow  through  the  sound 
field, 

e.  farfleld  radiation  pattern, 

f.  total  power  radiated  into  the  half-space. 

Nearfleld  holography  Is  based  on  conventional  holographic 
theory.  The  assumption  Is  that  a  group  of  localized  sources 
produces  a  wave  field  In  three  dimensions.  This  wave  field 
is  then  measured  across  a  two-dimensional  surface.  The 
amplitude  and  phase  information  recorded  in  two  dimensions 
Is  used  to  reconstruct  the  wave  field  throughout  the  entire 
half-space  located  above  the  sources.  This  can  be  done 
because  the  field  obeys  the  acoustic  wave  equation  and 
Green's  theorem  applies.  The  process  just  described  Is 
exact  in  theory,  but  In  practice,  conventional  holography 
poses  some  severe  restrictions  on  the  reconstructed  data. 
In  particular,  spatial  resolution  Is  usually  limited  by  the 


6 


radiated  wavelength.  In  other  words,  details  of  a  source 
located  less  than  a  wavelength  apart  cannot  be  resolved  and 
these  features  appear  smoothed  out.  This  makes  the  detailed 
study  of  sound  radiation  from  sources  operating  at  low 
frequencies  unfeasible  with  conventional  acoustical 
holography.  Conventional  techniques  are  also  limited  to 
reconstructing  a  scalar  'intensity'  proportional  to  the 
square  of  the  pressure.  It  Is  not  possible  to  determine  the 
vector  Intensity,  and  therefore  one  cannot  pinpoint  energy 
sources  or  construct  energy  flow  maps.  These  limitations 
are  due  to  experimental  restrictions  that  are  always  present 
In  optical  holography  but  do  not  necessarily  apply  to 
measurements  of  sound  waves.  Nearfleld  holography  removes 
these  experimental  limitations. 

The  sound  field  generated  by  any  monochromatic  source  can 
be  decomposed  Into  an  angular  spectrum  representing  plane 
waves  moving  In  different  directions.  Actually,  the  field 
Is  composed  of  propagating  plane  waves  and  evanescent  waves. 
Evanescent  waves  are  those  whose  amplitudes  decrease 
exponentially  with  distance  from  the  source.  This  decay  Is 
very  rapid,  with  amplitudes  typically  lOOdB  down  In  Just  one 
or  two  wavelengths  from  the  source.  This  Is  well  below  the 
dynamic  range  of  most  measuring  Instruments.  Conventional 
holographic  techniques  measure  the  wave  field  at  distances 
of  many  wavelengths  from  the  source  and  therefore  can  not 
record  the  evanescent  wave  Information.  The  propagating 
plane  waves  that  can  be  recorded  are  produced  by  features  of 
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the  source  that  vary  In  space  with  a  scale  of  length  that  Is 
greater  than  the  radiated  wavelength  and  this  results  In 
spatial  resolution  limitations.  Evanescent  waves  vary  In 
space  with  a  scale  of  length  that  Is  less  than  the  radiated 
wavelength  and.  If  measured,  allow  details  of  the  source  to 
be  resolved  regardless  of  the  size  of  the  radiated 
wavelength.  From  a  mathematical  point  of  view,  the 
evanescent  wave  Information  Is  needed  for  the  exact 
application  of  Green's  theorem.  An  example  of  the  dramatic 
effects  of  the  evanescent  wave  components  Is  Illustrated  In 
Flg.l  (Taken  from  Williams  and  Maynard,  Ref.  4)  which  shows 
various  aspects  of  a  wave  field  generated  by  two  simple, 
experimental  sources.  Fig. la  shows  the  pressure  field  In  a 
plane  very  close  (1/15  of  a  wavelength)  to  the  plane  of  the 
sources.  Even  this  close  there  Is  very  little  Information 
that  would  Indicate  two  distinct  sources.  Flg.Ib,  the 
reconstructed  pressure  using  conventional  holography,  shows 
the  effect  of  leaving  out  the  evanescent  wave  Information. 
The  details  of  the  two  sources  have  been  smoothed  over  a 
wavelength  and  any  Indication  as  to  the  location  and  number 
of  sources  Is  lost.  Flg.lc  shows  the  reconstruction  that 
retains  the  evanescent  information.  The  sources  can  be 
clearly  resolved. 

Nearfleld  holography  strives  to  measure  as  much  of  the 
evanescent  wave  Information  as  possible  by  using  an 
Instrument  with  a  large  dynamic  range,  and,  since  the 
evanescent  waves  decay  exponentially  with  distance  from  the 


Figure  1.  Holographic  Data  for  Two  Point  Sources,  (a)  Sound-Pressure 
Amplitude  in  the  Hologram  Plane,  (b)  Conventional  Holographic  Xe- 
tonstruction  of  the  Sources.  The  Resolution-Limiting  Wavelength  is 
Shown,  (c)  Reconstruction  of  the  Source  Intensity  Using  >;eariield 
Holography.  The  Two  Point  Sources  are  Clearly  Resolved. 
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source,  the  wave  field  should  be  measured  as  close  as 
possible  to  the  source.  Once  this  is  done,  the  evanescent 
as  well  as  propagating  wave  Information  can  be  recorded  and 
digitally  processed  with  an  exact  formulation  of  Green's 
theorem.  This  solution  seems  quite  elementary  and  the 
question  may  arise  as  to  why  this  was  not  done  previously. 
The  reason  seems  to  be  that  acoustical  holography  was 
developed  from  the  technology  of  optical  holography  where  It 
Is  impossible  to  measure  evanescent  waves  due  to  the  very 
small  optical  wavelengths.  In  low  frequency  sound  fields  It 
Is  easy  to  make  measurements  within  a  wavelength  from  the 
source  since  acoustic  wavelengths  can  be  a  meter  or  more, 
but  there  was  no  motivation  from  optical  holography  to  do 
so.  Additionally,  it  Is  not  obvious  that  measuring  the 
nearfleld  pressure  alone  Is  enough  to  characterize  how  a 
source  delivers  acoustic  energy  to  a  medium.  The  nearfleld 
pressure  (or  particle  velocity)  alone  Is  not  enough  to 
determine  how  a  vibrator  radiates  sound;  the  complete  vector 
intensity  Is  required  to  do  this.  If  nearfleld  measurements 
are  made,  an  exact  form  of  Green's  theorem  allows  one  to 
calculate  the  velocity  field  from  the  pressure  field  and 
then  obtain  the  vector  Intensity  from  these  two.  In  fact, 
any  desired  property  of  the  entire  sound  field  may  be 
calculated  with  this  technique.  Omitting  the  evanescent 
waves  precludes  exact  calculations  and  therefore  one  Is  left 
with  the  limits  of  conventional  holography.  Fortunately,  a 
large  class  of  radiators  can  be  measured  using  nearfleld 
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holography  and  a  great  deal  of  new  Infornatlon  becones 
available . 

1.2  Discussion  of  Theory 

An  overview  of  the  natheaatlcal  theory  used  In  nearfleld 
holography  Is  In  order  .  This  discussion  will  be  Halted  to 
single  frequency,  planar  sources  and  a  tlae  dependence  of 
exp(-lwt)  Is  assuaed.  The  fundaaental  equations  aake  use  of 
Green's  theorea  which  states  that  a  knowledge  of  the  coaplex 
pressure,  representing  aaplltude  and  phase,  across  an 
Infinite  plane  located  In  a  sound  field  allows  one  to 
propagate  to  any  other  parallel  plane  and  deteralne  exactly 
the  coaplex  pressure  there.  The  relationship  expressing 
this  Is  Rayleigh's  second  Integral: 

P(x,y,z)-C /2p<x' .y' I  2'«d}dx'dy'  (2) 

where  p(x',y',d)  -  the  pressure  In  a  source  plane 

located  at  z'  >  d 

p(x,y,z}  >  the  pressure  In  a  plane  z  away 
froa  the  source  plane 

3/3  z' [  exp(  lkR)/R]  I  ■  pressure-pressure  Green's 

function 

R  -  /  (x-x')2+(y-y')2  +  (2_2')2 
C  -  (-1/2’f  ) 

If  -  3.1416 
and  z  >  z' 


This  equation  allows  us  to  deteralne  the  pressure  at  a 
distance  z  froa  the  source.  But  we  would  like  to  be  able  to 
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■easure  the  pressure  In  a  plane  at  some  distance  z*h 
(hologram  plane)  and  use  the  mathematics  to  propagate  back 
to  the  source  plane  at  z>d.  To  do  this,  ve  rewrite  equation 
(2)  as: 

pCx.TiZ)  "  /2  P<x'.y'»d)  g(x-x' ,y-y' ,z-d)  dx'  dy'  (3) 

where  g(x-x" ,y-y' ,z-d)  -  C  3 / 3 z' [ exp( IkR) /R] | ^ 

This  form  Is ' recognizable  as  a  two-dimensional  convolution 
Integral  that  may  be  written  as: 

p(x,y,z)  -  p(x',y',d)  *  g(x-x' ,y-y' , z-d)  (4) 

where  *  denotes  the  convolution  operation.  The  convolution 
Integral  Itself  Is  quite  difficult,  but  by  employing  Fourier 
transform  techniques,  we  can  transform  the  equation  Into 
frequency  space.  Once  this  Is  done,  the  convolution  becomes 
a  simple  product  of  two  terms.  Once  the  product  Is  formed, 
the  Inverse  Fourier  transform  returns  the  relation  to  real- 
space  . 

The  two-dimensional  Fourier  transform  may  be  defined  as: 

P(kjj,ky,d)=  /2  p(x,y  ,d)exp(-lkjjx)exp(-ikyy)dxdy  (5) 

If  we  transform  both  sides  of  equation  (4)  we  have: 

■  P(lix,ky,<i)  G(k3j,ky,z-d)  (6) 

where  G  -  C  F{ 3/ 3z' [ exp( IkR) /R] | ^  } 


and  F  denotes  the  Fourier  transform 
operation. 
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Since  the  two-dimensional  transform  does  not  operate  on  z' , 
G  may  be  rewritten  as: 

G  -  C  3/3z'{Fte*p(lkR)/R] } I (7) 

With  some  work,  the  transform  of  the  Green's  function  Is 
found  to  be: 

expdz  ’^k^  -  k^2  _  V  2j 

F[exp(lkR)/R]  -  - - - ,  -/ -  (8) 

1  k^  -  kjj^  -  ky^ 

where  k“2Tr/X*2  ir  f/c 

k^  and  k  are  spa*‘lal  frequencies 
In  the  x^and  y  directions 

If  we  define  k^  •  k^-kjj^-ky^  ,  equation  (8)  becomes: 

F[exp(lkR)/R]  -  exp(lk2(z-d))/lk2 

After  taking  the  derivative,  the  transformed  Green's 
function  Is: 

G  -  exp(l  k^  (z-d)) 

and  the  final  form  of  equation  (6)  Is: 

P(k^,ky,z)  -  P(k^,ky,d)  expl Ik^(z-d) ]  (9) 

Equation  (9)  also  shows  the  effect  of  the  spatial 
frequencies  on  the  argument  of  the  exponential.  There  are 


two  cases: 
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1 .  k  ■  k^-k  ^-k  ^  ,  for  k^  >  k  ^  +  k  ^ 

*  *  ^  ^  y 

such  that  G  ■  exp[l(z-d)  k^-k^^-ky^  ] 

2.  k_  -  1  "k-^+k-^-k^  ,  for  k^  <  k-^  +  k„2 

such  that  G  ■  exp(-(2-d)  k_^+k_^-k^  ] 

*  7 

The  first  case  represeats  the  propagating  plane  waves  whose 
phases  change  with  distance  from  the  source.  The  second 
case  shows  an  exponentially  decaying  function  representing 
the  evanescent  waves.  The  function  decays  rapidly  with 
distance  and  suggests  that  aeasurement  of  the  evanescent 
waves  must  be  made  very  close  to  the  source.  Inverting 
equation  (9)  yields: 

P(k^,ky,d)  -  P(k^,ky,2)  expI-1  k^  (*“d)]  (10) 

Equation  (10)  allows  us  to  measure  the  pressure  In  some 
plane  above  a  source  and  'hack-propagate'  that  Information 
to  the  source  plane  located  at  z-d.  The  Inverse  Fourier 
transform  of  equation  (10)  returns  the  Information  to  real 
space : 

F"!  [  P(k^,ky,d)  ]  -  f’M  P(k3j.ky,z)  expt-1  k^  (z-d)]} 

p(x,y,d)  -  //  P(kjj,ky,z)  exp[-l  k^  (z-d)] 

*  exp(lk^x)  exp(lkyy)  dk^^dky  (11) 

Equation  (11)  Is  the  fundamental  relationship  of  nearfleld 
acoustical  holography.  It  should  be  noted  that  d  Is  not 
restricted  to  the  source  plane.  In  fact,  It  may  correspond 
to  any  plane  above  the  source  plane.  This  allows  us  to 
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determine  the  pressure  at  any  point  In  the  half-space 
located  above  a  source. 

The  vector  velocity  vhich  Is  necessary  for  calculating 
the  vector  Intensity  may  be  obtained  by  considering  that  the 
velocity  Is  proportional  to  the  gradient  of  the  pressure, 
1 .  e . 

^p  -  1  <JP  V  (12) 

where  ^  *  ( 2^^  3/  3x  +  Sy  ^  3y  +  8^  3/3  j] 

The  Fourier  transform  of  equation  (12)  yields: 

V(k^,ky,d)-(l/aiP)[kj^8^  +  kySy  "  l  ^  ^  3/  3  Z  ]  P(k^,ky,d) 

•  (  1/ up)  [kjjSj,  +  ^y®y  *^z*z  ^ 

*  exp[-l  k^  (z-d)]  (13) 

The  vector  velocity  Is  then  obtained  by  Inverse  Fourier 
transforming  equation  (13): 

v(x,y,d)  -  -C/up  [kjjSjj  +  kySy  +  k^S^j  P(k^,ky,z) 

*  exp[-l  k^  (z-d)] 

*  exp(lk^x)  exp(lkyy)  dk^  dky  (14) 

It  should  again  be  noted  that  z  Is  not  restricted  to  the 
source  plane  (z*d)  and  may  correspond  to  any  plane  above  the 

source.  The  modal  structure  of  the  vibrator  may  be  obtained 

by  considering  the  normal  component  of  the  velocity  at  the 
surface  of  the  source.  From  equation  (14),  the  normal 
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component  Is: 


v(x,y,<i)  -  -C/uP  /J  P(k^,ky,s)  expt-ik^(  z-d)  ] 

*  exp(ikjjX)  exp(lkyy>  dk,^  dky  (15) 

The  addition  of  a  time  factor  exp(-lvt)  allows  one  to  view 
the  motion  of  the  surface  through  time. 

We  have  derived  an  expression  for  the  pressure  In 
equation  (11)  and  the  vector  velocity  In  equation  (15),  and 
It  Is  now  possible  to  calculate  the  vector  Intensity,  I(r), 
from  equation  (1).  Since  the  wave  field  may  be 
reconstructed  In  any  plane,  mappings  of  the  Intensity  may  be 
produced  that  show  the  energy  flow  around  and  away  from  the 
source. 

Other  quantities  that  may  be  calculated  Include: 

The  power  radiated  Into  a  half-space, 

'if  ~  {1  Hx, 7, d)  dx  dy  (16) 

The  radiation  efficiency,  , 

-  W  /( p  c  S 

where  W  >  power  from  an  object  of  surface  area  S 

^  ■  average  RMS  velocity  of  radiating  surface 
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and  the  directivity  pattern, 

p(r,0,6)  -  i  uip  [exp(ikR)/R]  V(k,,ky,d)  (17) 

where  V  is  the  transformed  normal  velocity 
and  the  spatial  frequencies  are: 

k^  *  [2n/X]  cosQ  sin6 
ky  ■  [2ir/Xl  ain0  sln^ 

All  of  this  Information  may  be  obtained  from  one,  quick 
measurement  of  the  pressure  amplitude  and  phase  in  a  plane. 
The  time  savings  are  obvious  and  a  great  deal  of  information 
becomes  easily  available  to  the  researcher. 

1.3  Brief  Description  of  Measurement  Technique 

An  array  of  256  electret  microphones  are  suspended  in  a 
flat  plane  (the  hologram  plane  previously  mentioned)  and  are 
arranged  in  a  16  x  16  square  grid.  These  microphones  record 
the  sound  field  within  10  cm  of  the  source  plane.  The 
physical  spacing  between  the  microphones  may  be  varied  to 
change  the  resolution  of  the  system.  The  effective  sampling 
of  the  sound  field  is  enhanced  for  steady-state  sources  by 
scanning  the  radiated  sound  field  with  the  array  to  fill  in 
points  between  the  physical  spacing  of  the  microphones. 
This  is  actually  done  by  moving  the  source  on  a  two- 
dimensional  scanning  apparatus,  an  easier  operation  than 
moving  the  microphone  array.  The  source  is  usually  moved 
through  a  16  point  (4x4)  square  grid.  The  distance  between 
the  scanning  grid  points  is  determined  as  follows: 

d  ■  A/4  ,  where  d  «  distance  between  grid  points 

A  -  microphone  spacing  in  the  array 
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The  result  of  scanning  the  sound  field  Is  a  64  x  64  array  of 
data  points  rather  than  the  16  x  16  points  available  without 
scanning . 

The  microphones  are  sampled  to  give  a  time  sequence  for 
each.  Each  time  sequence  Is  Fourier  transformed  to 
determine  the  amplitude  and  relative  phase  at  each 
microphone.  These  values  are  stored  In  a  complex 
representation  and  used  with  the  processing  techniques 
described  earlier. 
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CHAPTER  TWO 


ARRAY  CONSTRUCTION 

2.1  Historical  Introduction 

When  this  project  began,  an  array  system  using  dynamic 
microphones  had  been  constructed  by  Williams  et  al^ .  This 
system  was  used  for  conventional  acoustical  holography 
measurements  of  steady  state  sources.  Shortly  after  I  began 
my  work.  It  was  discovered  Chat  Che  resolution  limit  of  a 
wavelength  that  Is  usually  associated  with  holography  need 
not  apply  In  some  acoustic  applications^.  This  new 
technique,  which  we  call  non-wavelength  limited  acoustical 

e 

holography  or  nearfleld  holography*',  required  that  Che 
measurement  be  made  quite  close  to  Che  source.  This  meant 
Chat  Che  measurement  plane  should  typically  be  within  10  cm 
(4  In)  of  Che  source  plane.  Much  work  went  Into  Che 
development  of  this  technique  ,  but  the  results  were  not  as 
good  as  we  originally  expected  they  would  be.  In  our  search 
for  an  explalnat Ion ,  we  realized  that  Che  problem  was  not  In 
the  theoretical  development,  but  rather  In  some  factors 
related  to  the  actual  physical  measurement.  First,  the 
sources  we  study  are  usually  electro-mechanlcally  driven. 
This  sets  up  a  varying  magnetic  field  around  the  source. 
The  dynamic,  or  moving  coll,  microphones  nearest  Che  source 
were  being  Influenced  by  this  field  and  therefore  were 
Inaccurately  recording  the  sound  field.  Second,  Che 


microphones  were  mounted  on  5  cm  high  by  .95  cm  wide  (2  In 
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by  3/8  la)  aluminum  beams.  The  refleccloa  of  sound  from 
these  beams  became  significant  when  the  source  was  moved 
close  to  the  array.  Finally,  the  microphones  were  not 
sensitive  to  the  lower  frequencies  that  the  new  technique 
would  allow  us  to  study.  All  of  these  factors  led  to  the 
decision  to  build  a  new  microphone  array  that  would  avoid 
these  problems. 

2.2  New  Array  Design  Considerations 

When  designing  the  new  array  system,  a  set  of  design 
goals  were  decided  upon.  These  are  discussed  In  the 
following  paragraphs. 

First,  the  microphones  used  in  the  new  array  should  not 
be  affected  by  nearby  magnetic  fields  and  have  a  broad,  flat 
frequency  response.  Inexpensive  electret  microphones  were 
chosen  because  they  fullfll  these  requirements. 

Second,  the  suspension  system  and  supporting  frame  should 
be  designed  to  securely  hold  the  microphones  In  position  and 
minimize  any  sound  scattering  and  reverberation  effects.  To 
this  end,  a  sturdy  aluminum  frame  ,  with  rows  and  columns  of 
thin  steel  piano  wires  strung  taughtly  across  It,  was 
designed.  Microphones  were  located  at  the  intersections  of 
the  wires. 

Third,  the  new  array  should  be  able  to  record  broad  band 


sources . 

To  do  this. 

a  sec  of 

time  samples 

from 

Che 

microphones 

Is  needed. 

Consider  a 

square  array 

of  N 

X  N 

elements . 

The  Ideal 

processing 
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be 

to 
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simultaneously  record  all  of  the  microphones  In  the  array 
analogous  to  taking  a  photograph  since  all  of  the 
microphones  would  be  recorded  at  the  same  Instant  In  time)  , 
but  this  would  require  Identical  processing  circuits  and 

would  be  quite  expensive.  Today's  state-of-the-art 
electronics  allow  processing  times  that  approximate 

simultaneous  recordings.  Therefore,  by  processing  a  single 
row  of  N  microphones  simultaneously  (with  N  Identical 
processing  circuits)  and  then  sequentially  processing  the  N 
different  rows  at  a  high  rate,  we  can  record  the  entire 
array  very  quickly.  Multiplexing  the  rows  with  analog 

switches  at  the  site  of  each  microphone  simplifies  the 
wiring  of  the  array  and  only  2N,  rather  than  N*” »  lengths  of 
cable  ace  required.  The  exposure  time  of  this  "acoustic 
photograph”  Is  determined  by  the  row  multiplexing  rate. 
Once  the  microphones  have  been  multiplexed,  analog  analysis 
of  the  sound  field  In  the  time  domain  Is  no  longer  possible. 
Therefore,  digital  conversion,  storage  and  analysis  Is 
necessary.  For  state-of-the-art  electronics,  the  data 
acquisition  rate  Is  limited  by  the  analog-to-dlgltal  (A/D) 
conversion  time;  a  typlcl  value  Is  l.Susec.  For  N*16  (a  256 
element  array)  the  minimum  exposure  time  of  such  an  acoustic 
photograph  Is  16  x  l.Susec  ■  24usec.  This  would  be  quite 
sufficient  for  studying  sound  sources  In  our  frequency  range 
of  Interest  as  will  be  discussed  later. 

The  signal  from  an  electret  microphone  Is  typically  in 
the  mllll-volt  region.  An  A/D  converter  requires  a  lOV 
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sigaal  for  proper  operation.  Therefore,  the  microphone 
signal  must  be  amplified.  If  the  amplifier  Is  used  after 
the  multiplexing  stage.  It  must  be  able  to  respond  at  the 
higher  multiplexing  frequency.  A  typical  operational 

amplifier  can  provide  60  dB  of  gain  at  the  higher 
multiplexing  rate.  If  a  larger  gain  Is  needed.  It  must  be 
provided  before  multiplexing  and  should  be  at  the  site  of 
the  microphone  to  eliminate  any  cable  noise.  A  calculation 
of  the  total  gain  needed  follows: 

Noise  level  for  a  typical  electret  microphone  -  0.2uV 
Signal  to  noise  ratio  for  good  reconstruction  >  30dB^ 
Signal  level  should  be  0.2  x  ■  6.3uV 

for  a  signal  level  of  lOV  at  A/B: 

Total  gain  -  20  Log[10/(6.3  x  10"^)]  *  124  dB 

Since  the  gain  after  multiplexing  can  be  60  dB, 
amplifiers  at  the  microphones  must  provide  64  dB. 

The  spacing  between  microphones  should  be  variable  to 
allow  different  aperture  sizes  and  resolving  capabilities. 
Using  N>16  rows  and  columns,  the  maximum  spacing  that  will 
fit  in  our  existing  quasl-anecholc  chamber  Is  17.8  cm  (7 
In).  The  minimum  spacing  that  allows  proper  mounting  of 
components  in  the  array  Is  2.54  cm  (1  in). 

To  determine  the  temporal  sampling  rate  and  highest 
recordable  frequency,  we  must  consider  the  Nyquist  sampling 
theorem  in  both  space  and  time.  For  a  microphone  spacing  of 
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R,  the  spatial  Nyqulst  theorea  allows  a  alnlaua  two- 
dimensional  wavelength  of  2R  to  be  accurately  recorded. 
Using  a  microphone  spacing  of  one  Inch,  the  highest 
allowable  frequency  Is  f  >  c/2R  *  6.6kHz.  If  any  higher 
frequencies  are  present,  spatial  aliasing  will  occur.  A  low 
pass  filter  must  be  used  before  multiplexing  to  avoid 
aliasing  In  the  time  domain.  The  minimum  time  between 
microphone  samples  Is  1/2(6600)  >  76usec,  so  the  entire  256 
element  array  must  be  sampled  within  this  time.  Thus  the 
obtainable  sampling  time  of  24usec  Is  more  than  adequate  for 
spatial  sampling  at  2.54  cm  Increments. 

The  digitized  data  flow  and  storage  of  the  microphone 
samples  must  now  be  considered.  The  flow  rate  Is  data 
words  In  76u8ec.  For  N  ■  16,  the  flow  rate  is  *  3.4  mega¬ 
words/sec,  too  fast  for  a  single  Input-output  (1/0)  channel. 
Therefore,  we  must  use  a  temporary  storage  unit  that  will 
accept  the  data  at  the  high  I/O  rate.  Since  there  are  N 
Identical  A/D  circuits  operating  In  parallel  for  the 
multiplexing,  let  each  circuit  have  Its  own  I/O  channel  and 
temporary  storage.  Then  the  flow  rate  becomes  N  data  words 
In  76usec,  or,  for  N  16,  200  kilo-words/ sec .  This  flow 
rate  is  possible  with  large  scale  Integrated  circuits. 

The  result  of  the  sampling  process  will  be  a  time 
sequence  for  each  microphone  In  the  array.  A  typical  length 
for  a  time  sequence  Is  1024  samples.  Taking  1024  acoustic 
photographs  with  the  N  x  N  array  would  result  In  a 
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'filmstrip'  of  duration  T  ■■  1024  x  76usec  ■  -  .08  sec  for 
our  6.6kHz  source.  Frequencies  below  this  allow  longer 
duration  filmstrips.  The  final  step  In  this  process  would 
be  to  read  the  data  (at  a  suitable  rate)  from  the  temporary 
storage  Into  permanent  storage  on  the  computer  for 
processing  by  reconstruction  programs. 

2.3  Interim  System 

Instead  of  constructing  16  A/D  and  temporary  storage 
circuits,  we  found  a  commercial  unit  that  performed  the  same 
operations.  It  utilizes  parallel  track-and-hold  circuits 
and  a  very  fast  A/D  converter  with  memory.  The  unit  Is  a 
LeCroy  CAMAC  Model  8210  Quad  10-8it  Transient  Digitizer  that 
can  do  four  simultaneous  A/D  conversions  at  a  rate  of  1  MHz 
and  store  the  Information  In  a  temporary  memory  unit.  The 
equipment  was  ordered,  but  before  It  arrived  we  constructed 
an  Interim  system  that  would  allow  us  to  continue  our 
nearfleld  studies  of  single  frequency  sound  sources.  The 
system  was  designed  to  work  with  the  LeCroy  unit  so  that  the 
transition  time  from  the  single  frequency  to  the  broad  band 
system  would  be  a  minimum. 

A  block  diagram  of  the  Interim  system  Is  shown  In  Figure 
2.  It  relies  on  an  AR-11  unit  In  our  PDP-Il/34  computer  to 
do  the  A/D  conversions  of  the  microphone  signals.  This  unit 
operates  at  a  rate  of  30  kHz.  Using  the  Hyqulst  theorem, 
the  highest  temporal  frequency  that  can  be  recorded  by  a  256 
element  array  with  a  30  kHz  sampling  rate  Is: 
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30,000  A/D  X  1  array  sweep  x  1  period  ■  58.59  periods 
sec  256  A/D  2  sweeps  sec 

-  58.59  Hz 

This  Is  much  Coo  low  to  be  useful.  To  study  higher 
frequency  sources,  we  decided  to  break  up  the  array  Into 
smaller  groups  Chat  we  call  blocks  and  take  our  acoustic 
photographs  one  block  at  a  time.  Processing  fewer 

microphones  at  one  time  raises  the  maximum  temporal 

frequency  Chat  we  can  record.  The  result  Is  a  time  sequence 
for  each  microphone  In  Che  array.  We  do  not  have  the 
capability  of  studying  broad  band  sources  since  all  of  the 

microphones  In  the  array  are  not  sampled  In  the 

'simultaneous'  manner  previously  described. 


Presently,  the  system 

Is  sec 

up  Co  continue 

our 

studies 

of  large 

(typically  30 

.5  X  61 

cm)  sources 

where  the 

wavelength 

of  the  radiated  sound 

is  typically 

larger  chan 

any  dimension  of  the 

source . 

We  chose 

the 

maximum 

microphone 

spacing  of 

17.8  cm 

Co  minimize 

any 

finite 

aperture  effects  with  Che  large  sources.  When  deciding  upon 
Che  low  pass  filter  breakpoint,  a  blockslze  of  nine 
microphones  (eight  array  microphones  and  one  reference 
measurement)  was  chosen  as  a  standard  size.  The  breakpoint 
calculation  Is  as  follows: 

30,000  A/D  X  1  block  meas.  x  1  period  ■  1667  Hz 
sec  9  A/D  2  block  meas. 


To  Insure  that  the  filter  was  active  at  this  value,  a 


ilok  Dicigram  of  liiturlm  Data  Ac(|iiiaLtion  Syutum. 
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breakpoint  of  1.5  kHz  was  selected.  The  length  of  the  time 
sequence  was  chosen  to  be  128  samples.  This  Is  a  good 
compromise  between  the  number  of  samples  required  for 
accurate  analysis  and  the  time  needed  to  sample  the  entire 
array . 

2.4  Details  of  Construction 


2.4.1  Microphone  Unit 

There  were  several  reasons  for  selec 
unit  rather  than  a  simple  microphone.  F 
voltage  from  an  electret  microphone  Is 
mllll-volt  region.  Amplification  of 

necessary  to  insure  good  transmission  thro 
wire  connecting  the  microphone  to 
electronics  with  a  minimum  noise  pickup, 
the  small  Field  Effect  Transistor  (FET)  In 
be  supplied.  Third,  locating  the  high  s 
microphone  allows  efficient  multiplexing 
reduction  In  the  required  length  of  conne 
previously  explained. 
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Mgiiru  3,  Ulock  Diagram  of  Microphone  Unit 


CHEMATIC  OF  MICROPHONE  UNIT  CIRCUIT 


Sciieinatic  of  Mlcro|>l>une  Unit  Circuit 


the  sound  field  Is  not  disturbed  significantly  by  the 
microphone.  The  microphone,  a  Prlmo  Co.  model  EM~60B, 
measures  10  mm  dia.  x  IS  mm  high  and  thus  satisfies  the 
small  size  requirement.  The  electret  microphones  were 
available  at  a  cost  of  72  cents  each.  Since  we  needed  256 
microphones  for  the  array  plus  some  spares,  this  price  was 
very  practical.  The  frequency  range  of  Interest  to  us  is 
from  100  to  2000  Hz  and  the  response  curve  for  these 
electrets  is  flat  through  this  region  as  well  as  above  and 
below  it.  The  FET  in  the  microphone  requires  a  -flOv  supply 
voltage.  Since  the  positive  supply  voltage  varies  between 
-M2  and  -flS  volts,  as  will  be  explained  subsequently,  a  lOv 
Zener  diode  is  used  to  maintain  the  proper  supply  voltage  to 
the  microphone. 

The  operational  amplifier  used  is  a  Texas  Instruments 
SN72558  dual  amplifier.  The  circuit  chosen  provides  30  dB 
of  gain  at  each  stage.  In  between  the  two  amplifiers  is  a 
swltchable  20  dB  attenuation  circuit.  This  is  used  if  a 
particularly  loud  source  is  being  studied  that  causes  the 
second  amplifier  to  clip.  Clipping  occurs  when  the  second 
amplifier  reaches  a  level  of  5  volts  peak  to  peak.  The 
schematic  shows  that  this  circuit  is  controlled  by  an 
electronic  switch.  The  switch, is  controlled  by  the  positive 
supply  voltage.  A  Zener  dio^e'  holds  the  gate  of  the  switch 
at  0  volts  when  the  supply  is  -MZ  volts.  This  corresponds 
to  the  switch  being  in  the  closed  position.  When  the  supply 
voltage  is  Increased  to  +15  volts,  the  gate  of  the  switch  is 
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supplied  with  4-3  volts  and  the  switch  opens,  activating  the 
attenuator  circuit. 

A  4~pole  filter,  2  poles  at  each  amplification  stage.  Is 
used  to  filter  out  any  high  frequency  components  that  are 
not  of  Interest  to  us.  This  prevents  any  aliasing  problems 
when  the  microphone  samples  are  Fourier  transformed.  The 
total  filter  rolloff  Is  24  dB/octave.  Presently,  the 
breakpoint  of  the  filter  Is  set  at  1.5  kHz.  Should  the  need 
ever  arise,  the  breakpoint  may  be  changed  by  a  plug-ln 
capacitor  board  on  the  microphone  unit 

The  high  speed  switch  Is  a  Harris  HI-200  SPST  CMOS  analog 
unit  that  contains  two  switches  In  one  package.  One  switch 
controls  the  attenuator  circuit  and  the  other  Is  used  to 
multiplex  the  microphones.  The  microphone  units  In  a  y 
direction  row  all  share  a  common  signal  out  line.  The 
switch  Is  used  to  turn  on  one  unit  at  a  time  along  this  y 
row . 

The  mounting  of  Che  components  Is  based  on  a  small 
printed  circuit  board  with  two  14-pln  DIP  sockets  attached 
CO  It.  This  can  be  seen  In  Figure  5.  A  socket  Is  located 
at  Che  lower  end  of  circuit  board.  This  Is  where  the 
microphone  plugs  Into  Che  unit.  The  operational  amplifier, 
filter  board,  and  switch  package  plug  Into  the  female  side 
of  Che  DIP  sockets.  The  top  six  pins  on  Che  male  side  of 
Che  DIP  sockets  are  used  to  plug  Che  unit  Into  sockets 
located  In  the  array  suspension  system.  These  six  pins 
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serve  Co  hold  Che  unlc  securely  In  place  and  make  Che 
necessary  elecCrlcal  connecClons  Co  Che  array  wiring.  If  a 
unlc  malfuncclons ,  Ic  Is  easily  replaced.  All  of  che  major 
componencs  on  Che  unlc  Icself  plug  In  and  are  also  easily 
replaced  If  chey  fall  Co  operaCe  properly.  As  wlch  Che 
microphone,  Che  microphone  unlc  muse  be  as  small  as 
possible.  To  Chls  end,  all  of  Che  componenCs,  capaclcors, 
reslsCors,  and  diodes,  were  chosen  by  chelr  size  and  were 
Che  smallesc  available  aC  Che  clme  of  consCructlon.  The 
final  measuremencs  of  a  compleCed  microphone  unlc  are  6  x 
1.5  X  1.5  cm . 

2.4.2  Microphone  Array  SupporC  Frame 

A  drawing  of  Che  microphone  array  supporc  frame  can  be 
seen  In  Figure  6.  IC  Is  consCrucced  of  10  x  7.5  cm  ( 4  x  3 
In)  aluminum  '1'  beam  and  measures  3.05  m  (10  fc)  square, 
measured  InCernally.  Two  of  Che  parallel  sides  spile  In  Che 
middle  so  chac  Che  frame  may  be  folded  In  half  and  moved. 
Two  3.25  m  (10  fC  8  In)  lengchs  and  four  1.625  m  (5  fC  4  In) 
lengchs  of  'I'  beam  were  used.  The  corners  were  formed  by 
making  a  45^  angle  cuC  aC  boch  ends  of  Che  long  pieces  and 
ac  one  end  each  of  Che  shore  pieces.  A  shore  lengCh  was 
welded  Co  each  end  of  Che  Cwo  long  lengCh  pieces.  The 
resulc  was  Cwo  'U'  shaped  pieces.  To  complece  che  frame, 
Che  cwo  'U'  pieces  were  Joined  cogecher  Co  make  a  square. 
The  spile  polnCs  are  held  cogecher  by  a  sandwich 
conf IguraClon  of  cwo  sceel  places  wlch  che  'I'  beam  In  che 
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middle.  The  plates  are  clamped  on  each  side  of  the  split 
pit  y  5  bolts.  This  securely  joins  the  two  frame  pieces. 
The  'I'  beam  size  was  chosen  because  of  Its  strength  and 
small  reflective  area.  The  frame  should  stay  as  square  as 
p  ssl  e  when  full  tension  Is  applied  to  the  piano  wires. 
The  worst  case  deflection  of  the  beam  was  calculated  as 
follows : 

Modulous  of  Elasticity  »  E  =10.3  x  10^  Ib/in^ 

Moment  of  Inertia  ■  I  «  1.31 

Load  (applied  at  center  of  beam)  >  U  >  16  x  301b  ■  4801b 

Length  of  Beam  «  1  >  120  In 

Max.  Oeflectlon  at  Center  ■  W  /  384  E  I 

“  0.16  In 

A  deflection  of  0.16  In  would  not  affect  the  frame  shape 
significantly. 

Holes  were  drilled  along  the  lower  Inside  edge  of  the  'I' 
beam  on  2.54  cm  (1  In)  centers.  These  holes  accept  tension 
adjusting  bolts  that  securely  hold  the  piano  wires.  The 
spacing  between  the  microphones  may  be  varied  from  17.8  cm 
(7  In)  to  a  minimum  of  2.54  cm  (1  in)  by  properly  locating 
the  tension  adjusting  bolts  In  the  holes.  This  allows 
adjustment  of  the  array-to-source  size.  Placing  the 
microphones  closer  together  Increases  the  resolution  of  the 
system. 

The  frame  Is  suspended  from  the  celling  by  a  rope  and 


pulley  system  attached  to  each  corner  of  the  array.  This 
allows  height  and  leveling  adjustment. 

The  piano  wires  need  to  be  securely  attached  to  the  frame 
and  stretched  so  that  they  do  not  sag  In  the  middle.  This 
was  accomplished  by  designing  tension  adjustment  bolts.  The 
bolts  must  be  long  enough  to  compensate  for  any  elongation 
In  the  wire  when  the  wire  Is  stretched.  The  following  Is  a 
calculation  of  the  wire  elongation  due  to  a  30  lb  tension 
force : 

F  =  Tension  Force  ■  30  lb 

L  ■  Wire  Length  *  120  In 

A  >•  Crossectlonal  Area  of  0.032  In  dla.  wire 
(0.032)2  /4  -  8.04x10"^  In^ 

E  *  Modulus  of  Elasticity  »  30x10^  Ib/ln^ 

Elongation  «  FL/AE 

■  0.15  In 

Thus,  the  change  In  length  Is  not  appreciable.  Two-Inch 
long,  3/8-lnch  diameter,  fully  threaded  bolts  were  selected. 
A  small  hole  that  would  accept  the  piano  wire  was  drilled 
perpendicular  to  the  bolt  head.  A  hole  for  a  set  screw  Chat 
firmly  clamps  Che  wire  In  Che  bolt  head  was  drilled  and 
Capped.  The  bolt  Is  placed  In  Che  'I'  beam  holes  and  held 
In  place  by  two  locking  nuts.  For  the  locking  nuts  to  work 
properly,  they  must  rest  against  a  parallel  face.  The 
Inside  edge  of  an  '1'  beam  Is  slightly  angled,  so  an  adaptor 
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p  "te  was  designed  that  provides  a  parallel  face  for  the  nut 
to  tighten  against.  It  also  serves  as  a  mounting  platform 
for  a  fuse  holder  and  wire  terminal.  The  complete  set  up 
can  be  seen  In  Figure  7.  A  worst  case  calculation  of  the 
sag  at  the  middle  of  the  piano  wire  was  made  to  determine  If 
the  suspension  system  would  hold  the  microphones  In  a  plane: 

0.11  Ib/mlke  unit  x  16  mike  units  ■  1.7  Ib^ 

New  Tensile  Force  In  Wire  -  {(30)^  +  (1.7)-^^ 

-  30.05  Ibf 

Elongation  Change  *  0.15  ln/30  Ib^  x  30.05  lb£ 

-  0.15025  In 

Length  Change  •  0.15025  -  0.15  •  2.5x10”^  in 

Deflection  Down  -  { (60.000125)^  -  (60)2>^'^^ 

-  0.12  in 

A  deflection  down  of  0.12  in  (calculated  for  the  worst  case 
of  all  16  microphones  being  placed  at  the  center  of  the 
piano  wire)  Is  not  appreciable  and  the  microphones  may  be 
considered  to  be  located  across  a  flat  plane. 

A  microphone  unit  plug-ln  socket,  an  eight  pin  DIP,  Is 
located  at  each  piano  wire  crossing.  Six  of  the  pins  are 
used  to  support  and  make  electrical  connections  to  the 
microphone  unit.  The  top  two  pins  were  removed.  Two  small 
holes  were  drilled  perpendicular  to  each  other  at  the  top  of 
the  socket.  The  piano  wires  thread  through  these  holes. 
When  the  piano  wires  are  tightly  stretched,  the  sockets  are 
oriented  In  the  same  direction.  This  assures  proper 
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allgnneat  and  placement  of  the  microphone  units  when  they 
are  plugged  In. 

2.4.3  Wiring  the  Array 

The  wire  used  to  make  electrical  connections  In  the  array 
must  be  as  small  as  possible  to  reduce  any  reflective 
effects.  To  this  end,  two-conductor  shielded  cable  with 
enamel-coated  28-guage  wire  was  chosen.  One  cable  runs 
along  each  piano'  wire.  In  both  directions,  making  three 
connections  to  each  socket  in  that  row.  The  positive  supply 
voltage  (+V)  and  signal  out  leads  run  In  the  y  direction  and 
the  negative  supply  voltage  (-V)  and  switch  gate  voltage  run 
In  the  X  direction.  It  Is  Important  that  the  gate  voltage 
and  the  signal  out  lines  do  not  run  In  the  same  cable 
because  of  possible  cross-talk  problems.  The  shield  In  both 
directions  Is  connected  to  the  ground  pin  of  the  plug-ln 
socket.  At  one  end  of  each  row,  a  supply  voltage  Is 
attached  through  a  protective  fuse  and  the  signal  out  or 
gate  voltage  Is  attached  to  a  mounted  terminal.  The  supply 
voltages  In  each  direction  are  bussed  together  and  attached 
to  the  Junction  box.  A  separate  line  Is  run  to  each  signal 
out  and  gate  voltage  row.  These  are  wired  to  two  ELKO 
multl-pln  connectors  that  plug  Into  a  Junction  box.  All 
cables  In  the  array  are  secured  with  wire  ties  to  keep  them 
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electronics  were  designed  by  Earl  Williams.  Any  one  of  the 
256  microphones  may  be  selected  for  output  by  an  eight-bit 
multiplexer,  the  HIKE  POSITION  unit  in  the  diagram.  As 
mentioned  earlier,  the  microphones  are  sampled  In  groups 
called  blocks.  Eight  switches  located  on  the  CONTROL  BOX 
select  the  size  of  the  block  as  Indicated  by  the  BLOCK  SIZE 
unit  of  Figure  8.  The  BLOCK  SIZE  number  does  not  change 
during  the  sampling  process.  When  the  POSITION  IN  BLOCK 
counter  receives  a  load  signal,  the  BLOCK  SIZE  is  latched 
Into  the  counter.  The  INITIAL  MIKE  POSTION  unit  stores  the 
position  of  the  lowest  numbered  microphone  In  the  block. 
When  the  hologram  data  taking  begins,  the  INITIAL  MIKE 
FOSIION  Is  set  at  zero,  which  Indicates  the  'first' 
microphone  In  the  array.  The  selected  MIKE  POSITION  Is  the 
sum  of  the  INITIAL  MIKE  POSITION  and  the  POSITION  IN  BLOCK 
counter.  Note  that  Initially  the  selected  microphone  Is  the 
sum  of  the  INITIAL  MIKE  POSITION  and  the  BLOCK  SIZE  and  that 
the  POSITION  IN  BLOCK  counter  Is  a  down  counter. 

At  the  beginning  of  the  data  taking  sequence,  the 
microphone  multiplexer  Is  disabled  and  a  reference  signal 
measurement  Is  taken.  Subsequently,  a  clock  pulse  enables 
the  multiplexer  and  the  first  microphone  Is  sampled.  The 


next  clock  pulse  triggers  the  POSITION  IN  BLOCK  counter  to 


41 


I 


It 
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count  down  by  one.  This  selects  the  next  microphone  In  the 
block  for  output.  The  countdown  process  continues  until  the 
POSITION  IN  BLOCK  Is  zero.  When  zero  Is  reached,  the  BLOCK 


SIZE  Is  again  latched 

Into 

the 

counter 

and 

the  data  taking 

Is  repeated.  The 

block 

Is 

sampled 

In 

this  manner  a 

selectable  number  of 

times 

to 

produce 

a  time  sequence  for 

each  microphone. 

When  the  units  discussed  above  are  configured  for  one 
particular  block  of  microphones,  the  INITIAL  MIKE  POSIION 
and  the  BLOCK  SIZE  are  summed  to  produce  the  NEXT  INITIAL 
MIKE  POSITION.  When  the  first  block  Is  finished,  the 
INITIAL  MIKE  POSITION  load  Is  activated  and  the  NEXT  INITIAL 
MIKE  POSITION  becomes  the  new  INITIAL  MIKE  POSITION.  This 
gives  a  new  first  microphone  position,  and  the  POSITION  IN 
BLOCK  counter  sequences  through  the  new  block. 

2.4.5  Details  of  Controlling  Electronics  Operation 
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Figure  9  Is  a  detailed  drawing  of  the  controlling 
electronics  and  Figure  10  shows  the  circuit  used  to  generate 
the  clock  pulses.  IC  19,  a  J-R  flip-flop,  enables  both  the 
count  down  system  and  the  multiplexer/reference  signal 
output.  q  starts  low.  This  Initiates  the  following: 

1.  signals  the  load  on  counters  IC  6  and  IC  8  to  latch 
the  BLOCKSIZE  from  the  FRONT  PANEL  SWITCHES  Into  the 


counters . 


Figure  9.  Detailed  Drawing  of  the  Data  Acquiaiticn 
Control  Electronics  Circuitry. 


uck  Circuit 
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disables  Che  taulclplexer  output  (through  ES)  and  closes 
switch  IC  13  to  place  the  REF  SIGNAL  IN  on  the  SELECTED  MIKE 
OUT  line. 

3.  disables  the  clock  from  counting  down  IC  8  through  the 
count-down  enable  NAND  gate  IC  15b. 

The  computer  signals  Che  clock  to  start.  The  first  clock 
pulse  causes  Che  following: 

1.  signals  the  computer  to  A/D  convert  the  signal  at  Its 
A/D  Input.  This  records  the  reference  signal. 

2.  toggles  IC  19  and  Q  goes  from  low  to  high  enabling  the 
count  down  through  IC  I5b.  IC  I6b  enables  the  multiplexer 
output  and  opens  switch  IC  13,  disabling  the  reference 
signal . 

The  microphone  selected  Is  determined  by  adders  IC  I  and 
IC  2.  They  sum  the  INITIAL  MIKE  POSITION  from  latches  IC  3 
a  and  b  with  the  counter  outputs  from  IC  6  and  IC  8.  IC  1 
signals  the  multiplexer  to  select  one  y-row  for  output.  IC 
2  signals  the  decoder  to  enable  a  gate  line.  This  selects 
the  first  microphone  In  the  block  for  output.  The  next 
clock  pulse  causes  the  following: 

1.  signals  the  computer  to  record  the  microphone  signal 

2.  decrements  counter  IC  8  by  one  placing  the  next 
microphone  In  the  block  on  line. 

Subsequent  clock  pulses  signal  counter  IC  8  to  count  down 
to  zero.  After  it  reaches  zero,  the  next  clock  pulse  sets 


all  outputs  high  and  BOR  goes  low. 


This  signals  IC  6  to 
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count  down  by  one.  The  last  microphone  in  the  block  Is 
selected  when  the  outputs  from  both  counters  are  all  low. 
The  next  clock  pulse  records  this  microphone  and  counts  down 
IC  8  from  all  low  to  all  high.  BOR  goes  low  and  IC  6  counts 
down  from  all  low  to  all  high.  The  BOR  on  IC  6  goes  low.  Is 
Inverted,  and  a  high  level  appears  at  pin  J  on  IC  19  causing 
Q  to  go  low.  The  first  microphone  In  the  block  Is  again 
latched  Into  the  counters  while  the  reference  signal  Is 
recorded  and  the  cycle  repeats,  sampling  the  block  a 
specified  number  of  times. 

When  the  sampling  of  the  first  block  Is  completed,  the 
computer  stops  the  clock.  Adders  IC  4  and  IC  5  sum  the 
BLOCK  SIZE  with  the  INITIAL  MIKE  POSITION  from  latches  IC  3 
a  and  b.  This  Is  the  NEXT  INITIAL  MIKE  POSITION.  The 
computer  signals  IC  3  t.>  Increment  the  block  and  the  NEXT 
INITIAL  MIKE  POSITION  Is  latched  to  the  output  of  IC  3.  The 
latch  output  Is  the  A  Input  of  adders  IC  1  and  IC  2  and  the 
new  block  Is  ready  to  be  sampled  by  the  counters. 
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CHAPTER  THREE 

ACCURACY  OF  ARRAY 

3 . 1  Purpo  se 

Once  the  construction  of  the  array  vas  completed  and  the 
microphone  data  were  being  properly  recorded,  It  was 
necessary  to  calibrate  the  microphones  and  determine  how 
accurately  the  system  was  reconstructing  the  sound  field 
being  generated  by  a  source.  We  wanted  to  know  the  limits 
of  the  system  and  consider  possible  sources  of  error.  This 
was  broken  down  Into  two  categories: 

1.  Reproducibility  -  determining  the  varlblllty  of  the 
reconstructions  between  successive  holograms. 

2.  Measurement  Accuracy  -*  compare  a  known  physical 
property  to  one  reconstructed  by  the  system  to  determine 
how  accurately  the  system  represents  that  quantity. 

Velocity  was  chosen  as  the  physical  quantity  because  the 
complex  velocity  at  the  surface  of  a  source  can  easily  be 
measured  by  an  accelerometer  and  directly  compared  to  the 
complex  velocity  holographically  reconstructed  In  the  source 
plane.  Verifying  the  velocity  rather  than  the  pressure  Is  a 
more  stringent  test  because  the  holographic  calculation  of 
the  velocity  Is  equivalent  to  taking  a  numerical  derivative 
of  the  data  which  greatly  magnifies  any  variations.  If  the 
velocity  Is  correctly  reconstructed,  then  the  pressure 
reconstruction  must  be  at  least  as  precise  since  the 


velocity  Is  obtained  from  the  pressure.  Knowing  that  the 
reconstructed  velocity  and  pressure  are  correct  insures  that 
the  reconstructed  intensity  is  correct  since  the  intensity 
is  the  product  of  the  in-phase  components  of  the  pressure 
and  the  velocity.  Therefore,  the  velocity  seems  to  be  a 
good  choice  in  determining  the  limits  of  the  system. 

3.2  Experimental  Technique 

The  microphone  calibration  procedure  will  be  reported  in 
detail  elsewhere  ,  but  a  brief  description  is  appropriate 
here.  A  calibration  chamber,  similar  to  a  piston  phone,  was 
designed  and  constructed.  A  reference  microphone  was  used 
to  determine  the  actual  complex  pressure  in  the  chamber  and 
all  of  the  microphones  were  calibrated  to  this  value.  The 
frequency  range  of  the  calibration  was  from  100  to  1500  Hz 
in  50  Hz  Increments.  This  proceedure  resulted  in  all  of  the 
microphones  being  calibrated  to  within  one  percent  amplitude 
and  one  degree  phase  of  each  other.  The  absolute  accuracy 
depended  on  the  accuracy  of  the  reference  microphone,  which 
was  estimated  to  be  about  three  percent. 

A  large,  flat  rectangular  plate  measuring  61  x  102  cm  was 
chosen  for  the  velocity  experiment  because  it  was  a 
relatively  simple  source  that  we  had  worked  with  previously. 
The  microphone  spacing  was  17.6  cm  yielding  a  4  by  4  square 
of  points  on  4.45  cm  centers.  A  grid  pattern  was  drawn  on 
the  plate,  symmetric  about  the  center,  that  corresponded  to 
points  in  the  reconstruction.  The  grid  consisted  of  lines 
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spaced  4.45  cm.  apart:  13  running  In  the  x-dlrection  and  25 
running  In  the  y-direction.  The  measurement  points  occurred 
at  each  x-y  crossing  giving  a  total  of  325  measurements. 

The  plate  was  mounted  on  the  scanning  apparatus  such  that 
the  grid  pattern  lined  up  as  closely  as  possible  with  the 
microphones.  The  plate  was  driven  at  248  Hz,  a  4,2  mode 
well  below  coincidence,  and  supported  by  foam  strips  along 
its  four  nodal  lines.  Two  holograms  were  recorded  in 
succession  with  care  taken  to  make  the  measurements  as 
similar  as  possible.  The  hologram  data  were  stored  in  the 
computer  for  analysis  at  a  later  time. 

The  plate  was  then  carefully  lowered  away  from  the  array 
making  sure  that  the  mounting  was  not  disturbed.  An  ENDEVCO 
Picomln  accelerometer.  Model  22,  was  used  to  record  the 
surface  motion  of  the  plate.  An  Ithaco  Dynatrac  3  Lock-In 
analyzer  was  used  to  determine  the  amplitude.  A,  and  the 
phase  angle,  8  ,  of  the  accelerometer  signal  relative  to  the 
place  drive  signal.  The  lock-in  output  was  two  DC  voltages 
proportional  to  A  sln6  and  A  cos6.  The  computer  recorded 
these  values  through  its  analog-to-digital  converter  and 
determined  Che  magnitude  and  phase  of  Che  plate  velocity. 
The  entire  plate  was  recorded  by  placing  Che  accelerometer 
at  each  of  Che  325  grid  points. 

Once  Che  reconstructed  and  measured  velocities  were 
scored  in  Che  computer,  they  were  compared  point  by  point  to 
determine  Che  error  between  the  two  velocity  files.  A 


49 

computer  program  was  written  to  plot  histograms  of  this 
error  In  the  following  way.  The  program  finds  the  maximum 
velocity  In  the  reconstruction  and  sets  up  20  velocity  bins. 
Each  bln  la  1/20  th  of  the  maximum  velocity  In  width.  The 
computer  then  scans  through  the  reconstruction  file  and 
finds  all ~bf  the  velocity  values  that  fall  within  each  bln. 
For  each  value  In  a  bln,  the  error  between  the  reconstructed 
and  measured  velocities  Is  calculated.  The  averaged  error 
of  all  the  velocity  values  that  fall  within  a  bln 
constitutes  the  total  error  for  that  bln.  The  error  Is  then 
plotted  as  a  function  of  the  20  velocity  bins. 

An  Error  Index,  a  number  that  would  represent  the  overall 
error  of  the  histogram,  was  determined  as  follows: 


Error  Index  * 


SUMrEi*Wi*m 

SUM[Wi*Nl] 


where:  El  ■  the  average  percent  error  In  bln  '1' 

Wi  *  the  weighting  factor  for  bln  '1' 

N1  >  the  number  of  values  In  bln  '1' 


The  weighting  factor  takes  the  following  form: 


Wl  ■  number  of  bins  -1+1;  1<  1<20 


where  1  «  highest  velocity  bln 
20  >  lowest  velocity  bln 


Because  we  expect  the  higher  velocity  regions  to  be  more 
Important  In  the  sound  field  reconstructions,  their  errors 
are  assigned  a  higher  weighting  factor. 

To  eliminate  any  offset  error,  the  velocity  values  were 


normalized  by  taking  the  ratio  of  the  reconstructed  to 
measured  values  at  the  center  of  the  plate  and  multiplying 


so 

the  measured  velocities  by  this  ratio.  This  greatly 
Improved  the  file  comparison  and  lowered  the  Error  Index. 
To  determine  If  this  was  the  best  normalization  procedure, 
an  Iteration  technique  was  used.  The  normalization  was 
Incremented  from  some  starting  value  In  small  steps  until 
the  lowest  Error  Index  was  obtained.  This  normalization  was 
then  used  for  the  histogram  plot. 

Once  the  histogram  data  taking  and  analysis  had  been 


completed,  we 

decided  to  run 

a  quick  check 

of  the  array 

calibration  to 

determine  If 

the  system 

was  correctly 

recording  the 

pressure  field. 

A  reference 

microphone  was 

placed  In  the 

array  plane. 

Measurements  of 

the  pressure 

over  a  point  on 

a  source  were 

recorded  by  the 

reference  and 

an  array  microphone.  The  errors  between  these  measurements 
were  compared  at  different  frequencies  and  were  found  to  be 
Inconsistent.  We  found  that  the  discrete  Fourier  transform 
(DFT)  was  returning  a  pressure  amplitude  that  was  always 
lower  than  the  reference  amplitude.  This  error  was  caused 
by  the  fact  that  the  source  was  being  driven  at  a  frequency 
F^  that  did  not  correspond  to  any  of  the  discrete 
frequencies  (multiples  of  Z’^/Ndt,  where  the  quantities  N  and 
dt  are  described  below)  set  up  by  the  DFT.  It  Is  easy  to 
show  that  to  correct  this  problem,  one  must  multiply  the 
amplitudes  returned  by  the  DFT  by  a  correction  factor  of  the 


f  orm : 
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Correction  Factor  ■  SIN  (x) 

(x) 

where  x  ■  Tr*N*dt*(Fjj  -  F^) 

N  ~  number  of  samples 
dt  ■  sample  time 

Fjj  ■  nearest  DFT  bln  frequency  to  source  frequency 
F^  ■  source  frequency 

This  factor  was  not  originally  Included  In  the  programs 
because  we  were  not  Interested  In  the  absolute  accuracy  of 
the  system  at  that  time.  Incorporating  this  correction 
factor  Into  the  data  taking  programs  solved  the  varying 
error  problem,  resulted  In  a  typical  error  between  reference 
and  system  measurements  of  about  four  percent  or  0.3  dB,  and 
Improved  the  absolute  accuracy  of  the  system.  The  relative 
accuracy  remained  the  same. 

3.3  Results 

The  reproducibility  of  the  system  was  determined  by 
taking  two  successive  holograms  and  using  the  histogram 
program  described  earlier  to  compare  them.  The  Error  Index 
from  this  histogram,  seen  In  Fig. 11,  shows  that  the  system 
Is  typically  reproducible  to  within  seven  percent  or  about 
0.59  dB. 

A  histogram  used  to  determine  the  measurement  accuracy 
can  be  seen  In  Fig. 12.  The  Error  Index  for  this  histogram 
shows  that  the  system  error  Is  typically  within  eighteen 
percent  or  1.^  dB.  This  would  seem  to  Indicate  that  the 
system  Is  correctly  representing  the  velocity  over  the 
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Magnitude  Error  in  Percent 


Showing  Measurement  Accuracy 
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surface  of  the  source,  but  sone  of  the  absolute  values  of 
the  velocity  differ  from  those  measured  by  the 
accelerometer.  The  phase  differences  between  the  hologram 
and  measured  data  are  typically  less  than  5  degrees. 

Other  possible  sources  of  error  that  time  did  not  allow 
me  to  study  are  as  follows: 

1.  The  standing  wave  pattern  set  up  within  the 
measurement  space  changes  when  the  source  Is  moved 
beneath  the  array  to  fill  In  points  betwen  the  actual 
microphone  spaclngs.  Changing  the  standing  wave  pattern 
may  change  the  pressure  recorded  by  the  microphones  and 
therefore  Introduce  an  error  Into  the  data.  The  effects 
of  changing  the  standing  wave  pattern  would  be  quite 
difficult  to  determine,  but  the  problem  may  be  eliminated 
by  keeping  the  source  stationary  and  moving  the  array  to 
fill  In  points.  Standing  waves  could  still  exist  within 
the  space,  but  the  pattern  would  not  change  appreciably 
If  the  array  was  moved.  To  do  this,  a  system  for  moving 
the  array  would  have  to  be  designed  and  constructed. 

2.  The  hologram  data  are  smoothed  during  the 
reconstruction  process  to  make  the  reconstruction  less 
noisy.  The  smoothing  process  changes  the  amplitude  and 
phase  of  the  data  slightly,  but  Is  necessary  for  clear, 
analyzable,  reconstructions.  This  must  also  be  a  source 


of  error  and  requires  a  detailed  study  to  fully 


55 


understand  Its  effects* 

3.  The  microphones  In  Che  array  are  not  all  exactly  on 
17.8  cm  centers.  This  was  realized  after  these 
experiments  were  completed.  Some  microphones  may  be  as 
much  as  one-half  Inch  offset  from  their  neighbor.  The 
result  of  this  error  Is  that  hologram  points  may  not  be 
lining  up  exactly  with  accelerometer  data  points  on  Che 
plate . 

4.  The  positioning  system  may  not  be  moving  the  source 
orthogonally  Co  the  array,  again  resulting  In  Che 
hologram  points  not  lining  up  with  accelerometer  data 
points.  If  Che  array  was  moved  Instead  of  the  source, 
this  source  of  error  would  no  longer  exist. 


CHAPTER  POUR 


EXAMPLE  OP  TECHNIQUE  USING  A  RELATIVELY  SIMPLE  SOURCE 
4.1  Introduction 

The  rectangular  plate  used  for  the  verification 
experiments  discussed  In  Chapter  3  nicely  demonstrates  the 
technique  In  actual  usage.  The  center  driven  plate  Is  a 
relatively  simple  source  In  that  there  are  no  ribs  or 
discontinuities  across  It  and  Its  modal  patterns  are  well 
known.  Therefore,  we  can  predict  vliat  the  reconstructed 
characteristics  should  be  and  compare  them  with  the  actual 
reconstructions. 

When  setting  up  the  experiment  for  Chapter  3,  a  mounting 
for  the  plate  had  to  be  selected.  This  required  that  the 
plate  be  supported  In  the  same  manner  throughout  the 
repeated  application  and  removal  of  the  accelerometer.  The 
easiest  mounting  scheme  that  would  fulfill  this  requirement 
was  to  place  the  plate  on  a  one- Inch- thick  sheet  of  foam. 
This  would  uniformly  support  the  plate  and  not  allow  it  to 
move  during  the  experiment.  To  Isolate  the  foam  from  the 
scanning  apparatus,  two  two-inch- t hick  sheets  of  fiberglass 
were  placed  between  them.  Both  the  foam  and  the  fiberglass 
were  cut  to  the  dimensions  of  the  plate.  Reconstructions  of 
this  configuration  showed  that  there  was  a  significant 
velocity  developed  in  the  air  Just  outside  the  boundaries  of 
the  plate,  which  was  probably  caused  by  diffractive  effects. 


I 
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The  data  smoothing  process  tends  to  average  the  velocities 
recorded  at  the  plate  edge  with  the  velocities  recorded  in 
the  air  just  outside  the  edge.  Since  the  air  velocity  is 
usually  180  degrees  out  of  phase  with  the  plate  edge,  this 
effect  would  pull  down  the  velocity  amplitude  measured  at 
the  edge.  It  should  be  noted  that  this  averaging  would 
happen  even  in  the  theoretical  limit  of  zero  velocity 
outside  the  plate  boundaries.  We  wanted  to  minimize  this 
effect  and  decided  to  try  a  different  mounting. 

Knowing  that  the  plate  was  being  driven  in  a  (4,2)  mode, 
we  decided  to  remove  the  foam  sheet  and  replace  it  with  foam 
strips  located  along  the  four  nodal  lines  across  the  short 
dimension  of  the  plate.  This  configuration  supported  the 
plate  correctly  and  we  hoped  that  this  would  reduce  the 
diffractive  effects.  Reconstructions  showed  that  the 
velocity  outside  the  plate  had  diminished  but  was  still 
present . 

Next,  we  thought  that  the  large  reflecting  area  of  the 
fiberglass  sheet  directly  below  the  plate  might  be 
contributing  to  the  diffractive  effects.  We  decided  to 
remove  these  sheets  and  replace  them  with  support  strips 
consisting  of  4  x  4  inch  blocks  of  fiberglass  extending 
along  the  short  plate  dimension.  The  foam  strips  used 
before  were  placed  on  cop  of  Che  fiberglass  and  were  again 
located  along  Che  nodal  lines.  This  configuration  supported 


the  plate  correctly  and  further  decreased  the  diffractive 
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effects.  It  was  therefore  selected  as  a  satisfactory 
mounting.  Previous  work  with  the  plate  point-supported 
Indicated  that  the  diffractive  effects  were  less,  but  the 
point-supported  plate  would  not  have  been  stable  enough  for 
the  accelerometer  measurements. 

A. 2  Results 

Figure  13  shows  an  overhead  view  of  the  plate  velocity 
with  the  full  f Iberglass/foam  sheet  mounting.  The  symbols 
In  the  legend  represent  levels  that  are  within  5,  10,  15, 
and  20  dB  less  than  the  maximum  level  recorded  on  the  plate. 
The  positive  (POS)  and  negative  (NEC)  symbol  sets  refer  to 
areas  of  vibration  that  have  a  phase  difference  of  180 


degrees  In  the 

velocity 

reconstruc 

tlons 

or  to 

areas  of 

positive  and 

negative 

Intensity 

In 

the 

Intensity 

reconstructions. 

By  comparing  this 

f  Igure 

with 

Figure  14 

(which  shows  the  plate  velocity  with  support  strip 
mounting) ,  one  can  see  that  aside  from  a  lessening  of  the 
diffractive  effects,  the  nodal  structure  and  velocity 
amplitude  are  quite  similar.  However,  when  comparing  the 
Intensities  of  these  two  support  cases,  an  Interesting 
effect  Is  observed.  The  Intensity  reconstruction  of  the 
foam  mounted  plate  Is  shown  In  Figure  15.  Note  the  large 
positive  Intensity  region  at  the  center  of  the  plate  where 
the  driver  Is  located.  Figure  16  shows  the  Intensity 
reconstruction  for  the  strip  supported  plate.  Comparing 
these  two  reconstructions  reveals  a  large  change  In  the 
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Intensity.  The  strip  supported  plate  shows  a  large  negative 
Intensity  region  near  the  driver  with  neighboring  areas  of 
positive  Intensity.  A  discussion  of  this  effect  follows: 

The  real  part  of  the  Intensity  Is  defined  as  one-half  the 
product  of  the  In-phase  components  of  the  pressure  and  the 
velocity.  1  .e. ; 

I  -  1/2  p  V  coa(  0  ) 

where  9  ^  the  difference  In  phase  angles 

between  the  pressure  and  the  velocity 

Analysis  of  the  foam-supported  plate  data  reveals  that  6 
near  the  center  of  the  plate  Is  typically  -68  degrees.  The 
cosine  of  -68  Is  positive  and  therefore  a  positive  Intensity 
Is  observed.  The  strip  supported  plate  shows  a  typical  9  of 
-98  around  the  center  with  neighboring  areas  having  a  9  of 
about  -81.  This  small  change  In  9  results  In  a  large 
Intensity  change.  A  similar  change  In  9  for  the  foam 
supported  plate  would  not  result  In  as  great  an  Intensity 
change.  Since  9  for  the  strip-supported  plate  Is  close  to 
-90,  the  sound  field  Is  quite  reactive.  This  almost  totally 
reactive  sound  field  does  not  produce  much  radiated  energy 
and  Is  therefore  a  poor  radiator  of  sound.  The  sound  field 
around  the  foam  supported  plate  Is  much  less  reactive, 
making  It  a  better  radiator  of  sound  even  though  the  plate 
Is  subjected  to  much  more  damping.  This  may  be  explained  by 
considering  the  radiation  from  an  Infinite  plate  that  Is 
point  driven.  The  flexural  waves  produced  by  the  driver 
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propagate  away  from  the  driven  area  through  the  plate. 
Since  the  plate  la  Infinite,  these  flexural  waves  are  not 
reflected.  In  this  case,  the  driven  area  produces  the 

O 

radiated  energy”.  In  our  case,  the  foam  supporting  the 
finite  plate  can  be  expected  to  absorb  some  of  this  energy. 
If  enough  of  the  flexural  wave's  energy  Is  absorbed  by  the 
foam,  the  damped  plate  begins  to  act  as  the  Infinite  plate 
and  most  of  the  radiated  energy  Is  produced  at  the  driver. 

On  the  other  hand,  once  most  of  this  damping  Is  removed 
by  strip-supporting  the  plate,  the  flexural  wave's  energy  Is 
not  as  highly  attenuated  and  more  of  this  energy  may  be 
reflected  from  the  plate's  boundaries.  This  results  In  a 
more  reactive  field  near  the  plate  and  less  radiated  energy 
Is  produced.  There  are  still  some  losses  present  and  some 
sound  energy  Is  radiated  as  can  be  seen  In  Figure  17.  This 
figure  nicely  shows  the  non-radlatlng ,  circulating  energy 
flow  present  In  the  strip-supported  plate. 
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CHAPTER  FIVE 

EXAMPLE  OF  TECHNIQUE  USING  A  COMPLEX  SOURCE 
5.1  la croductloa 

In  the  previous  chapter,  results  of  a  plate  source  with 
different  mountings  were  presented.  The  technique  was  able 
to  detect  changes  In  the  Intensity  due  to  different 
mountings  where  there  seemed  to  be  no  corresponding  change 
In  the  velocity.  The  next  and  final  step  In  this  work  was 
to  use  the  technique  to  study  a  complex  radiator.  To 
qualify  as  being  complex,  a  vibrator  may  have  a  non-trlvlal 
shape  or  Inhomogeneous  composition,  and  an  exact 
mathematical  solution  for  Its  sound-generating  properties 
does  not  exist.  To  enhance  the  relevance  of  this  research, 
It  would  be  helpful  If  the  source  had  some  general  academic 
Interest.  To  this  end,  we  chose  a  string-excited  enclosed 
cavity  with  radiating  port,  better  known  as  a  guitar,  to  be 
the  complex  source. 

The  guitar  Is  Indeed  a  complex  vibrator.  Its  shape  Is 
not  readily  adaptable  to  any  coordinate  system  and  the 
Internal  ribbing  makes  the  plates  Inhomogeneous.  The  source 
has  been  of  Interest  In  the  past^'  and  our  technique 
seems  Ideally  suited  to  study  It  since  we  can  reconstruct 
many  Important  properties  of  the  source  with  one 
measurement.  Previous  work  gives  us  an  Idea  of  what  to 
expect  from  the  plate  velocity,  but  to  the  author's 
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knowledge  there  have  never  been  any  studies  undertaken  to 
experimentally  determine  the  acoustic  energy  produced  and 
radiated  by  the  guitar. 

S.2  Experimental  Technique 

A  Martin  steel-stringed  acoustic  guitar  was  chosen  for 
the  experiments.  Unlike  many  other  experiments,  the  guitar 
was  totally  Intact  and  tuned  to  concert  pitch.  The  guitar 
was  point  supported  at  three  locations:  the  neck  and  two 
points  on  the  backside  edge  of  the  lower  bout  (The  guitar 
body  Is  referred  to  as  the  bout.  The  lower  bout  Is  the  area 
below  the  rose,  or  sound  hole,  and  the  upper  bout  Is  the 
area  above  the  rose.  In  this  discussion,  the  bout  will 
refer  to  the  top  of  the  guitar  unless  otherwise  noted).  In 
an  effort  to  excite  the  Instrument  as  closely  as  possible  to 
actual  playing  conditions,  the  steel  strings  were  used  to 
drive  the  Instrument.  A  small  electro-magnet  was  placed 
near  the  string  to  be  excited.  An  AC  current  was  passed 
through  the  electro-magnet  at  the  frequency  necessary  to 
oscillate  the  string  at  resonance.  An  electro-mechanical 
shaker  placed  on  the  bridge  of  the  guitar  would  excite  plate 
modes  more  dramatically  than  the  driven  string,  but  we  were 
Interested  In  examining  the  guitar  as  It  would  be  played. 


The  actual  study  was  divided  Into  two  parts.  The  goal  of 


the  first  study  was  to  locate  modes  of  the  guitar  top  plate 
and  compare  various  results  obtained  with  our  technique  to 


those 


techniques  (holographic 


obtained 


with 


other 
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Interferometery  la  particular).  The  second  study  focussed 
on  the  behavior  of  the  guitar  In  between  these  so  called 
natural  modes.  By  obtaining  reconstructions  of  the  top 
plate  over  the  first  two  octaves  of  the  Instrument's 
frequency  range,  we  hoped  to  see  changes  In  the  top  plate 
velocity  and  radiation  characteristics.  This  Information 
has  not  previously  been  available  and  should  be  of  Interest. 

5.3  Results 

5.3.1  First  Study 

At  first,  we  wanted  to  be  able  to  reproduce  some  of  the 

I 

results  found  In  the  literature.  There  appear  to  be  several 
accepted  modes  of  vibration  for  the  top  plate  of  the  guitar. 
These  can  be  seen  In  Figure  18.  The  first  plate  mode 
appears  as  a  breathing  mode  where  the  entire  lower  bout  Is 
vibrating  In  phase.  The  second  plate  mode  appears  as  a 
dipole  mode  where  the  lower  bout  breaks  up  Into  two  areas 
vibrating  180^  out  of  phase.  I  was  able  to  locate  both 
modes  In  the  experimental  guitar  but  will  concentrate  on  the 
dipole  mode  In  this  discussion.  The  frequency  of  the  dipole 
mode  was  determined  by  placing  a  shaker  on  the  bridge  and 
two  accelerometers  near  points  on  the  plate  that  could  be 
expected  to  be  vibrating  180^  out  of  phase.  Once  the 
frequency  of  the  dipole  mode  had  been  determined,  a  string 
was  tuned  to  that  frequency  and  excited  by  the  electro¬ 
magnet.  Velocity  reconstructions  at  this  frequency  reveal 
the  strong  dipole  mode  In  the  lower  bout,  but  other  effects 
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are  present  as  can  be  seen  In  Figure  19  (The  symbols  In  the 
legend  represent  levels  that  are  within  5,  10,  IS,  and  20  dB 
less  than  the  maximum  recorded  on  the  top  plate.  The 
positive  (POS)  and  negative  (NEC)  symbol  sets  refer  to  areas 
of  vibration  that  have  a  phase  difference  of  180  degrees  In 
the  velocity  reconstructions  or  to  areas  of  positive  and 
negative  Intensity  In  the  Intensity  reconstructions).  Here 
the  motion  of  the  air  In  the  rose  Is  evident  as  well  as  some 
diffractive  effects  outside  the  boundaries  of  the  guitar. 
These  diffractive  effects  are  probably  due  to  the  motion  of 
the  back  plate  of  the  guitar.  It  Is  Important  to  note  that 
the  Interferometry  technique  used  to  create  the  pictures  In 
Figure  18  cannot  represent  the  motion  of  the  air  In  the  rose 
or  outside  the  guitar  boundaries.  A  dipole  source  Is 
generally  accepted  to  be  a  poor  radiator  of  sound.  Using 
the  results  obtained  from  Interferometry,  one  might  be  led 
to  conclude  that  the  guitar  does  not  radiate  well  at  this 
frequency.  However,  our  reconstructions  show  an  area  of 
large  velocity  located  at  the  rose.  The  Intensity 
reconstruction  shown  In  Figure  20  Indicates  that  the  rose  Is 
a  substantial  energy-producing  area.  The  energy  flow  map  of 
Figure  21  clearly  shows  that  the  rose  Is  the  major 
contributor  to  the  radiated  sound  field.  Figure  22  shows 
the  nearfleld  short-circuit  that  occurs  due  to  the 
circulating  energy  flow  at  t  <wer  bout  of  the  guitar. 
This  results  in  little  radiated  energy  originating  at  the 
lower  bout.  Being  able  to  record  the  motion  of  the  air  In 
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X  (inches) 

KJKiir*i  20.  IiiLuiisilty  Rucoastructloii  of  CitJCdr  Dipole  Mode. 

Note  That  Most  of  the  Energy  la  Produced  at  the  Rose. 


WAVELENGTHS 

l'’lgure  21.  Eiturgy  t’low  Mai>  of  (Uiitar  Dipole  Mode  Taken  Across  Center  of  the  l.ong  Dimension. 
Note  I'low  Into  I’age  Over  the  Bridge. 


l'‘lfeiire  22.  blueir^y  l-'luw  Map  of  CiiJtar  Dipole  Taken  Across  Lower  Bout 

Tliroiitili  Bridge.  Note  That  Circulating  flow  Creates  a  Near- 
fleld  Short-Circuit. 
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the  rose  becomes  very  important  in  the  study  of  how  a  guitar 
produces  sound  energy. 

5.3.2  Second  Study 

The  preliminary  work  indicated  that  the  plate  modes  of 
this  guitar  were  reasonably  similar  to  published  results  on 
other  guitars.  Additionally,  more  information  is  available 
with  this  technique  in  that  we  are  able  to  record  the  motion 
of  the  air  as  well  as  that  of  the  plate.  The  next  step  was 
to  study  the  guitar  over  a  range  of  frequencies  to  observe 
its  behavior.  The  frequencies  chosen  were  chose 

corresponding  to  notes  available  in  a  tuned  guitar.  Each 
string  was  excited  at  Che  open,  fret  2,  and  fret  4 
positions.  Eighteen  data  files  resulted  covering  a 

frequency  range  of  79  to  425  Hz,  or  Che  notes  E2  to 

sharp.  The  strings  were  clamped  at  frets  2  and  4  by  a  capo 
(a  common  guitarist's  Cool  for  raising  Che  base  tuning  of 
Che  guitar).  The  electro-magnet  was  used  to  excite  Che 
selected  string,  and  Che  other  strings  were  damped  to 
eliminate  any  sympathetic  excitation. 

The  velocity  reconstructions  are  shown  in  Figures  23  and 
24.  The  cavity  resonance  for  this  guitar  is  around  100  Hz. 
At  Che  cavity  resonance,  one  would  expect  the  guitar  Co  be 
acting  as  a  Helmholtz  resonator.  The  air  in  the  rose  would 
be  acting  as  a  piston,  and  therefore  we  would  expect  to  see 
this  effect  in  our  reconstructions.  This  is  indeed  Che  case 
as  shown  in  Figure  23.  But  examination  of  the  other 
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frequencies  show  chat  Che  air  In  Che  rose  Is  always  an  area 
of  large  velocity  over  this  range.  This  effect  Is  alssed  by 
techniques  such  as  Interferometry.  In  fact,  even  at 
frequencies  with  large  plate  velocities,  as  at  160  Hz,  Che 
rose  Is  still  quite  dominant. 

The  first  plate  resonance  appears  around  170  Hz  and  this 
Is  shown  clearly  In  Che  178  Hz  case  of  Figure  24.  But  the 
plate  resonance  does  not  span  a  large  frequency  Interval  and 
at  189  Hz  Che  rose  Is  again  dominant.  This  appears  to  be 
true  up  to  298  Hz,  where  the  lower  bout  Cakes  over  and 
becomes  more  significant  Chan  the  rose  velocity.  Figures  23 
and  24  seem  to  Indicate  that  there  Is  more  Chan  one 
frequency  Chat  excites  the  breaching  mode  In  an  assembled, 
string  excited  guitar.  The  final  3  frequencies  of  Figure  24 
nicely  show  the  evolution  of  the  dipole  mode  from  a 
breathing  mode  at  400  Hz.  It  Is  also  Interesting  to  note 
the  motion  of  Che  neck  at  certain  frequencies  In  Figure  24. 
Subsequent  accelerometer  measurements  confirmed  Chat  this 
was  a  real  effect.  This  small  area  certainly  does  not 
contribute  to  the  radiated  sound  field  directly  but  may  h&ve 
some  effect  on  Che  width  of  resonance  peaks  due  Co  Its 
coupling  to  Che  body  of  the  guitar. 

Figure  25  Is  a  plot  of  the  ratio  of  Che  volume  velocities 
from  Che  rose  and  the  Cop  place  over  Che  frequency  range 
studied.  The  dotted  line  indicates  a  ratio  equal  to  I, 
corresponding  to  equal  volume  flow  from  Che  rose  and  Che  top 
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plate.  A  ratio  greater  than  one  Indicates  that  the  rose 
dominates  the  volume  flow  and  a  ratio  less  than  one 
Indicates  that  the  top  plate  dominates  the  volume  flow  (The 
volume  flow  Is  calculated  by  summing  the  velocity  over  that 
region).  This  plot  shows  that  below  300  Hz,  the  rose  Is 
responsible  for  most  of  the  volume  flow  produced  by  the 
guitar.  Between  300  and  400  Hz  the  top  plate  becomes  more 
significant,  and  at  425  Hz  (the  dipole  mode  previously 
discussed)  the  rose  Is  again  dominant.  The  frequencies 
above  400  Hz  are  about  two  octaves  above  the  cavity 
resonance  and  this  may  explain  the  Increase  In  activity  at 
the  rose.  A  study  of  frequencies  above  425  Hz  Is  necessary 
to  determine  whether  or  not  the  rose  continues  to  be 
significant  or  If  the  top  plate  again  dominates. 

The  Intensity  reconstructions  for  this  study  are  shown  In 
Figures  26  and  27.  As  with  the  velocity,  the  rose  Is 
dominant  up  to  266  Hz.  At  189  Hz,  the  lower  bout  shows  an 
area  of  negative  Intensity,  but  apparently  the  rose 
generates  enough  energy  to  produce  a  significant  net 
radiation.  If  this  negative  area  could  be  eliminated,  the 
guitar  would  surely  radiate  better  at  this  frequency. 
However,  changes  In  the  guitar  that  would  eliminate  the 
negative  area  at  189  Hz  may  produce  problems  at  other 
frequencies.  This  Is  a  complex  problem  In  need  of  further 
s  tudy . 

The  Intensity  reconstruction  at  364  Hz  is  interesting 


tensity  Reconetriicti 
the  frequency  Range 
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because  the  rose  seems  to  be  much  less  significant.  The 
lower  bout,  while  it  has  a  small  negative  area,  seems  to  be 
the  energy-producing  area  at  this  frequency.  Quite 
possibly,  the  rose  and  this  negative  area  are  acting  as  a 
nearfleld  short-circuit  as  was  discussed  in  the  first  study 
(Time  did  not  allow  verification  of  this  effect  through 
energy  flow  mappings).  If  this  were  the  case,  the  lower 
bout  would  be  acting  as  the  major  energy  producing  area. 
This  is  a  curious  effect  since  at  all  other  frequencies,  the 
rose  is  a  dominant  area.  A  close  study  of  frequencies  near 
364  Hz  may  reveal  a  frequency  where  the  rose  is 
substantially  reduced  in  both  velocity  and  intensity.  The 
negative  area  seen  at  364  Hz  may  expand  to  give  another 
dipole  mode.  Evidence  to  support  this  was  obtained  in  the 
preliminary  study  where  a  dipole  mode  without  a  large 
contribution  from  the  rose  was  observed.  This  may  be  seen 
in  Figures  28  and  29. 

A  plot  of  the  ratio  of  the  power  from  the  rose  to  the 
power  from  the  plate  is  shown  in  Figure  30.  This  plot  again 
reflects  the  dominance  of  the  rose  as  an  energy-producing 
area  below  300  Hz.  The  plate  dominates  between  300  and  400 
Hz  and  the  rose  is  again  the  dominant  area  at  425  Hz. 
Care.ul  examination  of  the  volume  velocity  and  power  ratio 
plots  of  Figure  25  and  Figure  30  show  that  they  are  quite 
similar.  The  intensity  seems  to  follow  the  volume  flow  in 
this  guitar  at  frequencies  below  425  Hz.  This  is 
significant  because  the  volume  flow  is  easier  to  measure 
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than  the  Intensity.  Accelerometers  can  determine  the  plate 
velocity,  and  pressure  measurements  at  two  closely  spaced 
points  perplndlcular  to  the  plane  of  the  top  plate  may 
determine  the  velocity  and  resulting  volume  flow  in  the 
rose.  If  this  were  true  for  other  guitars.  It  could  be  an 
easier  method  of  finding  the  power  produced  by  a  guitar  over 
this  range  of  frequencies. 

5.4  Discussion  and  Conclusions 

Finally,  the  results  from  this  study  should  be  compared 
to  results  published  by  other  authors.  This  author  feels 
that  using  an  Intact,  string-excited  guitar  is  a  more 
realistic  approach  to  studying  the  Instrument,  but  published 
results  obtained  with  such  a  guitar  have  not  been  found.  An 
understanding  of  the  natural  modes  of  the  plates  Is 
Important,  but  analyzing  the  behavior  of  the  assembled 
system  Is  necessary  to  better  understand  the  Instrument. 
The  work  that  I  would  like  to  concentrate  on  Is  that  of 
Caldersmlth^ ^  and  Firth^^. 

The  first  point  that  I  would  like  to  address  Is  the 

notion  that  below  the  cavity  resonance,  the  guitar  is  a  poor 

radiator  because  the  motion  of  the  air  in  the  rose  is  180® 

1 4 

out  of  phase  with  the  motion  of  the  top  plate.  Firth 
claims  that  at  frequencies  below  the  cavity  resonance,  the 
rose  does  not  enhance  the  radiation  of  sound  from  the  guitar 
and  that  the  radiated  SPL  falls  off  at  a  rate  of  18 
dB/octave.  An  analogous  acoustical  circuit  Is  presented  to 
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than  the  intensity.  Accelerometers  can  determine  the  plate 
velocity.  and  pressure  measurements  at  two  closely  spaced 
points  perplndlcular  to  the  plane  of  the  top  plate  may 
determine  the  velocity  and  resulting  volume  flow  In  the 
rose.  If  this  were  true  for  other  guitars.  It  could  be  an 
easier  method  of  finding  the  power  produced  by  a  guitar  over 
this  range  of  frequencies. 


5.4  Discussion  and  Conclusions 


Finally,  the  results  from  this  study  should  be  compared 
to  results  published  by  other  authors.  This  author  feels 
that  using  an  Intact,  strlng~exclted  guitar  Is  a  more 
realistic  approach  to  studying  the  Instrument,  but  published 
results  obtained  with  such  a  guitar  have  not  been  found.  An 
understanding  of  the  natural  modes  of  the  plates  Is 
Important,  but  analyzing  the  behavior  of  the  assembled 
system  Is  necessary  to  better  understand  the  Instrument. 
The  work  that  I  would  like  to  concentrate  on  Is  that  of 
Caldersmlth^ ^  and  Firth^^. 


The  first  point  that  I  would  like  to  address  Is  the 
notion  that  below  the  cavity  resonance,  the  guitar  Is  a  poor 
radiator  because  the  motion  of  the  air  In  the  rose  Is  180^ 
out  of  phase  with  the  motion  of  the  top  plate.  Flrth^^ 
claims  that  at  frequencies  below  the  cavity  resonance,  the 
rose  does  not  enhance  the  radiation  of  sound  from  the  guitar 
and  that  the  radiated  SPL  falls  off  at  a  rate  of  18 
dB/octave.  An  analogous  acoustical  circuit  Is  presented  to 
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support  this  claim.  Caldersmlth^ ^  claims  that  at 
frequencies  below  the  cavity  resonance,  air  flows  freely  In 
and  out  of  the  sound  hole  In  antiphase  with  the  top  plate 
motion  Implying  a  dipole  effect  between  the  rose  and  the  top 
plate  In  support  of  Firth's  claim. 

Caldersmlth  models  the  guitar  as  a  reflex  enclosure.  The 
body  and  back  are  represented  as  rigid  elements  with  a 
flexable  plate  and  sound  hole  on  top.  At  low  frequencies, 
this  system  Is  mass  controlled  so  an  Inward  motion  of  the 
plate  compresses  the  air  In  the  cavity  and  results  In  an 
outward  flow  at  the  rose.  Caldersmlth  Justifies  the  use  of 
this  model  by  measuring  the  velocity  of  the  top  and  back 
plates  of  a  guitar  at  the  top  plate  fundamental  resonance 
frequency.  It  Is  not  surprising  that  the  top  plate  motion 
at  this  frequency  could  be  greater  than  the  back.  The 
conclusion  drawn  Is  that  the  motion  of  the  back  plate  of  the 
guitar  may  be  Ignored  when  studying  the  lower  range 
performance  of  the  guitar. 

During  the  course  of  the  second  study,  accelerometer 
measurements  of  the  top  and  back  plates  were  recorded  for 
some  frequencies  to  compare  with  the  above  reported  results. 
One  such  data  set,  taken  at  79  Hz,  revealed  that  the  upper 
and  lower  bout  of  the  top  plate  were  vibrating  180^  out  of 
phase  resulting  in  a  low  net  volume  flow.  But  the  velocity 
of  the  air  In  the  rose  was  about  10  times  greater  than  that 
of  the  top  plate  at  this  frequency.  How  then  Is  this  large 
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volume  flow  at  the  rose  produced?  Accelerometer 

measurements  across  the  back  plate  revealed  that  the  back 
was  vibrating  In  phase  and  had  an  amplitude  2  to  3  times 
greater  that  the  top  at  this  frequency.  This  motion  was 
found  to  be  In  phase  with  the  air  motion  In  the  rose. 
Obviously  this  Is  where  the  large  volume  flow  at  the  rose  is 
produced.  This  Indicates  that  the  motion  of  the  back  plate 
Is  significant  and  should  not  be  Ignored.  Velocity 
reconstructions  at  frequencies  above  and  below  the  cavity 
resonance  show  little  motion  In  the  top  plate.  Therefore, 
we  must  conclude  in  contrast  to  Caldersmlth  and  Firth,  that 
the  back  plate  must  be  a  very  Important  mechanism  In  the 
production  of  sound  at  these  lower  frequencies  and  should  be 
considered  In  subsequent  studies. 

The  second  point  I  shall  address  is  that  of  top  plate 
modes  In  the  guitar.  It  appears  that  the  modes  generally 
accepted  to  be  natural  modes  of  the  top  plate  may  occur  at 
several  frequencies  In  a  string  excited  guitar.  When  a 
plate  is  shaker  driven,  one  would  expect  a  mode  that  Is 
characteristic  to  that  plate  to  be  excited  at  only  one 
frequency.  This  is  the  observed  trend  In  shaker  drlva 
guitar  plates^^*  The  results  presented  In  this  study  cf 

a  string  excited  guitar  Indicate  that  these  modes,  or.  what 
may  be  slight  variations  of  these  modes,  occur  at  more  than 
one  frequency.  This  leads  to  the  conclusion  that  a  complex 
relationship  between  the  strings,  ribs,  body,  and  neck  of 
the  guitar  causes  similar  modes  at  different  frequencies. 
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Adding  these  elenents  (strings,  ribs,  etc.)  inezsases  the 
number  of  degrees  of  freedom  and  turns  single  resonances 
into  multiplets  with  normal  modes  similar  to  those  of  the 
simple  plate.  One  would  think  that  a  more  or  less  flat 
response  curve  with  many  overlapping  multiplets  rather  than 
a  few  sharp  resonances  would  be  desirable  because  this  would 
eliminate  any  notes  that  sound  much  louder  than  their 
neighbor . 

The  next  logical  step  in  this  research  would  be  to  study 
the  back  plate  of  the  guitar  in  a  manner  similar  to  that 
used  in  the  top  plate  study.  This  should  lead  to  a  better 
understanding  of  the  interaction  between  the  two  plates. 
The  development  of  an  analogous  acoustical  circuit 
containing  the  back  plate  would  certainly  assist  in  this 
endeavor  and  should  be  realizable  after  studying  the  back 
plate.  The  effects  of  the  neck  being  coupled  to  the  body 
and  its  resultant  vibrations  are  also  of  interest.  Of 
course  the  ribbing  structure  is  very  important.  There  has 
been  some  work^^  in  which  the  vibrations  of  the  top  plate  of 
guitars  with  different  numbers  of  ribs  were  examined.  It 
would  be  interesting  to  study  any  changes  in  the  intensity 
with  the  change  in  the  number  of  ribs.  Lastly,  the  effects 
of  the  sides  of  the  guitar  should  be  investigated.  It  seems 
that  the  sides  are  generally  assumed  to  be  rigid,  but  when 
the  guitar  is  played,  one  can  feel  the  sides  vibrating. 

A  useful  research  team  would  consist  of  a  luthier  and  a 
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scientist  using  nearfleld  holography  as  a  research  tool* 
The  luthier's  knowledge  of  guitar  construction  would  allow 
him  to  make  subtle  changes  In  the  guitar's  construction. 
The  effect  of  these  changes  could  be  studied  by  nearfleld 
holography,  and  a  better  understanding  of  how  structural 
changes  affect  radiation  could  be  achieved.  Taking  a 
particular  Instrument  that  Is  considered  to  have  superior 
qualities  and  determining  Its  radiation  characteristics  over 
a  frequency  range  of  Interest  could  create  a  'standard  of 
excellence'.  The  luthler/sclentlst  team  could  endeavor  to 
create  an  Instrument  that  matched  the  radiation 
characteristics  of  this  standard  of  excellence  and  should 
thereby  develop  an  Instrument  of  exceptional  quality.  In 
mass-produced  Instruments,  nearfleld  holography  could  be 
used  as  a  quality  control  tool  for  Instruments  In  higher 
price  ranges.  At  any  rate,  the  holographic  technique 
presented  here  could  act  as  a  major  contributor  to  the 
understanding  of  musical  Instrument  radiation  and  possibly 
raise  the  standard  of  quality  In  many  applications. 
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